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We show that the recent observation of the time modulation of two-body weak decays of heavy 
ions reveals the mass content of the electron neutrinos via interference patterns in the recoiling 
ion wave function. From the modulation period we derive the difference of the square masses 
Am 2 pa 22.5 x 10 -5 eV 2 , which is about 2.8 times larger than that derived from a combined analysis 
of KamLAND and solar neutrino oscillation experiments. It is, however, compatible with a data 
regime to which the KamLAND analysis attributes a smaller probability. 

The experimental results imply that the neutrino mixing matrix violates unitarity by about 10%. 

PACS numbers: 



I. INTRODUCTION 

At the GSI in Darmstadt, the experimental storage ring ESR permits observing completely 
ionized heavy atoms I or hydrogen-like heavy ions Ih over a long time [3, and thus to measure 
the time dependence of their weak two-body decays Ih — * I + Ve or I — ► Ih + Ve- The first is the 
well-known electron-capture (EC) process. The virtue of such experiments is that the properties 
of the neutrino or antineutrino can be deduced from measurements of the time dependence of the 
transition observing only the initial and final ions. The special efficiency of these experiments 
becomes clear in the Dirac sea interpretation of the second process, where the initial ion simply 
absorbs a negative-energy antineutrino in the vacuum. Since the vacuum has all negative-energy 
states filled, the vacuum is a source of negative-energy neutrinos of maximally possible current 
density, i.e., the best possible neutrino source in the universe. This is why the ESR experiments 
yield information on neutrino properties with great precision even if the targets and exposure 
times are quite small, in particular much smaller than the 2.44 x 10 32 proton-yrs (2881 ton-yrs) 
in the famous KamLAND experiments [H, which are only sensitive to the much less abundant 
positive-energy neutrinos produced by nuclear reactors. 

Apart from the neutrino mass difference, the experiment reveals also another important property 
of the presently popular neutrino mixing scheme: the matrix which expresses the neutrino flavor 
states into fixed-mass states must be nonunitary to explain the data. The measurement determines 
the degree of nonunitarity to be roughly 10%. 

II. TWO-NEUTRINO MIXING 

To illustrate this we consider here at first only the two lightest neutrinos. According to Pon- 
tecorvo 0, IH, the Dirac fields of the physical electron and muon-neutrinos Vf = {y e ,v^), the 
so-called flavor fields , are superpositions of neutrino fields i/j = (v\, v-i) of masses m\ and mi'- 

v e (x) = v\{x) cos + i>2{x) sin#, v^x) =— Vi{x) sin 9 + v% (x) cos9 , (1) 

where 9 is a mixing angle. This is, of course, the neutrino analog of the famous Cabibbo mixing 
of up and down quarks. The free Dirac action has the form 

A = J2 I d 4 xP f (x) (i'fd^-M)uf(x), 
f J 

where 7^ are the Dirac matrices, and M. is a mass matrix, whose diagonal and off-diagonal elements 
are nif = (m^m^) and m e ^ = m^ e , respectively. The eigenvalues m, = (mi, 7712) are related to 
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TO/ 



by ilia, 



toi cos 2 9 + mi sin 2 8, m^ = mi sin 2 9 + mi cos 2 i 



m efl = TO Me = (to 2 — to | ) sin cos 0. (2) 
The weak transition between the electron e and its neutrino v e is governed by the interaction 

Ant =5 J d 4 xW-{x)J + ^(x)+h.c. =g J d i xW;{x)e{x) 1 ^{\- lb ) Ve {x)+h.c, (3) 

where 75 is the product of Dirac matrices £7777. 

Since the interaction ([3|) involve only the flavor fields ([1]) , the states of masses to^ will always be 
produced as coherent superpositions. The weakness of the interaction will allow us to calculate the 
shape of the mixed wave packet from perturbation theory. Consider the decay / — > In + v e which 
is a superposition of the states of masses m\ and mi. The formulas will be applicable for electron 
capture if we exchange Mh by the mass M of the bare ion and deal with outgoing neutrinos. 

In the center-of-mass (CM) frame of the initial bare ion of mass M, the final H-\ike ion has 
the same momentum as the antineutrino Vi (i = 1,2), whose energy is uJi = coki.i — \/k 2 + m 2 
determined by 



M = M H + Q = uo l + ^M 2 H + k 2 =u Jl + ^M 2 H +uo 2 -m 2 1 z = 1,2, (4) 

so that 

uj, = [(2M H + Q)Q + m 2 ]/2{M H + Q). (5) 

Subtracting loi and u)% from each other we find the energy difference 

ml - m\ Am 2 

Alj = W9 - wi = = . 

2 1 2M 2M v ' 

The denominator M is of the order of 100 GeV and much larger than Ato 2 , so that Aui is extremely 
small. It is the difference of the recoil energies transferred to the outcoming ion by the antineu- 
trinos of masses toi and mi. Without recoil, we would have found the four orders of magnitude 
larger energy difference at the same momentum Awt = Wk,2 — ^k,i = (Am 2 +w^ 1 ) 1 '' 2 — u>u.,i ~ 
Am 2 /2ajk,i ~ Ato 2 /2Q. This is the frequency with which the incoming negative-energy neutrino 
current of momentum k oscillates in the vacuum. Note that although Auj is small, the momentum 
difference Ak = ki — k\ associated with the energies oj\^i is as large as Auk, but has the opposite 
sign. 



III. EXPERIMENTS 

The best experimental results are available for the EC-processes reported in Ref. [lj, where an 
electron is captured from the K-shell and converted into an electron-neutrino which runs off to 
infinity. On the average, the decay is exponential with a rate expected from a standard-model 
calculation. In addition, however, the decay rate shows modulations with a frequency Aw. The 
experimental results are [1] 

i4o Pr 58+ ^ i4o Ce 58+ . A w « 0.890(11) sec^ 1 (Q = 3 386keV), (7) 
'^Pm 6 ^ -> ^ 2 Nd 60+ : Alu » 0.885(31) sec" 1 (Q = 4470keV). (8) 

In both cases, the period of modulation is roughly 7 sec, and scales with M (see Fig. 1). The decay 
rate has the form A(i) = A(0)(1 + acos(Aut + Ac/))) with a modulation amplitude of a — 0.18(3). 

We expect theses modulations to be associated with the frequency Aw of Eq. ©, and thus to 
give information on Ato 2 . Inserting the experimental numbers for Alu into Eq. ([6]) and taking into 
account that the particles in the storage ring run around with 71% of the light velocity with a 
Lorentz factor 7 1.43, we find for both processes @ 



Ato 2 « 22.5 x 10~ 5 eV 2 . 



(9) 
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Figure 1: Modulations of decay rate for the processes 14 5 gPr 5S+ -> "gCe 584 " and 1 gPm 60+ -> 1 ^Nd 60+ . 
The period is in both cases roughly 7 sec. The inserts show the frequency analyses. Plots are from Ref. [J]. 
The decay rate is modulated by a factor 1 + a cos(Aujt + A0) with a — 0.18(3). 
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Figure 2: The upper KamLAND regime of 2006 [lOJ is compatible with our result Am 2 « 22.5 x 10" 5 eV 2 . 



This is by a factor « 2.8 larger than the result Am 2 » 7-58±g;|x 10 _5 eV 2 favored by the KamLAND 
experiment [1, [ll| , but it lies close to their less favored result , which the authors excluded by 
2.2a in 2005, and now by 6a (see Fig. 2). 

So far we do not yet understand the origin of this discrepancy. One explanation has been 
attempted in Ref. [12| where the authors investigate the influence of the strong Coulomb field 
around the ion upon the process. 



IV. ENTANGLED WAVEFUNCTION 



For a theoretical explanation of the modulations, we first simplify the situation and ignore all 
spins and the finite size of the ions. Then the decay of the initial ion I into the ion Ih plus an 
electron-antineutrino v e . can be described by an effective interaction for this process is 



d 4 x I H (x)v e {x)I(x) =, 



cos 



f5/| r (ar)i/i(a;)7(x)+sinf97^(x)i^(a;)/(a;)] , (10) 



where I(x), v e {x), and Ih(x) are the field operators of the involved particles. In the CM frame, 
the initial ion is at rest, the final moves nonrelativistically. The outgoing wave is spherical. The 
role of the antineutrino creation operators in v\ and vi is simply to create a coherent superposition 
of two such waves with the two different k- and w-values calculated above. The combined outgoing 
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wave function will be 

/ /^[cos0e-^P t+ ^|/H(p))e-^ lt+4kx |^(k)) 

+ sm6»e-' iB P* +lpx |/ jFf (p))e~ l " k ' 2 * +lkx |z>2(k))] . (11) 

The states |Pi(k)), |z/2(k)), in turn, can be reexpressed in terms of the electron- and muon-neutrino 
states as 

|Pi(k)) = cos#|P e (k)) - sin0|^(k)), |j/ 2 (k)) = sin0|P e (k)) + cos#|^(k)). (12) 
Thus we find for the transition to an electron-neutrino of any momentum the effective action 

d 3 k f 

— ^(Pef-kJI^tatlO) = gj d 3 xI H {x)I{x)v^{x) (13) 
with a spacetime-dependent potential 

/d 3 k 
[cos 2 O e *"*.i*+a~ + sin 2 9e^ t+ *x] . (14) 

In Born approximation we find from this the scattering state of the recoiling ion Ih- 

MrpW-t)^ = ~ 9 - fcos 2 9e ,(fcir - Wlt » + sin 2 6 e^-^A . (15) 

This wave carries a radial current density of ions Ih 
a 2 

ir° = ITT — o [cos 4 6ki + sin 4 9 k 2 + sin 2 9 cos 2 9 (ki + k 2 ) cos(Afc r - Auj t)} . (16) 
Mht 

In order to find the decay rate we integrate this over a sphere of radius R surrounding the initial 
ion, choosing for R any size <C 1/A/c ss 10 4 m. For this surface we find the outgoing probability 
current density 

k 

P = 4tt 5 2 — [l - i sin 2 (29) + \ sin 2 (20) cos(Aw t)\ , (17) 

where we have approximated k\ and ki by their average k. 

This P can explain directly the observed modulations of the decay rate of the initial ions. 
The is only one problem: the amplitude of modulations are predicted to be a — ^ sin 2 (2#)/[l — 
■jsin 2 (2#)] = 0.72. Experimentally, however, a is much smaller. It has the value 0.18(3). 

The discepancy is explained by a missing contribution to the deay rate. So far we have only 
included the contribution of the effective action (fTOf . There is, however, also a second effective 
action which is generated by the matrix elements 

d 3 k ( 

(P M (-k)|^ int |0) =g / d 3 xI H {x)I{x)v^(x) (18) 



(27T) 

where 

(2ir) 

This can be derived directly from Eqs. (fl"2"|) and (fT3")) . Here the Born approximation yields the 
scattering state of the recoiling ion Ih ■ 

(x|V> (+) ;t)^ = -sin (9 cos 6> h***'—"-*) - e *(fe»r-wat)l _ ( 20 ) 
r L 

Its radial current density is now 

2 

3r" = TJ—2 sin2 6 cos2 + fe 2) [1 - cos(Afc r-Aw t)] . (21) 

Now we have another problem: the modulations of this current cancel the modulations of the 
current (fTT|) . We may suspect that this has to do with the fact that there are three neutrinos 
which we must take into consideration. 
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V. THREE-NEUTRINO MIXING 



Let us now include all three known neutrinos v ei v^,v T . Their fields are denoted by v a with 
a = e, /J,, t. These fields are combinations of three fields with definite mass 

v a = U a iVi, cr = (e,/i, t), (22) 

The mixing matrix U„i is called Maki-Nakagawa-Sakata matrix, or short MNS-matrix, the neutrino 
analog of the 3x3 Cabibbo-Kobayashi-Maskawa matrix for the mixing of the quarks d, s, b, It is 
commonly assumed to be unitary, i.e., to satisfy the relation 

Y,KjU„i=h- (23) 

Its standard parametrization is the following product of four simple unitary matrices 

10 \ / ci3 s 13 e- lS \ ( ci2 sia \ / e iai 
U = ( c 23 s 23 1 - 8t2 c 12 e ia * I (24) 

-s 23 c 23 / V -*i 3 e l5 0ci3/\0 1 J \ 01 

where = sinfljj, = cos#jj. For quarks, the unitarity relation (I23p is presently in the focus of 
experimental and theoretical studies in many research groups 14]. For leptons, the data have so 
far been insufficient to test it. 

Generalizing (fTB"]) . (fl~4f and (fl8l) . (fT9|) . we have from each flavor cr an effective potential 

H 3 k f 

' (9 a (-k)\A int \0)=g d 3 xI H (x)I(x)v p Jx) (25) 



(27T) 



with a potential 

v*. (*) = J (03 E ^^•^ mk *. (26) 

This produces an outgoing wave function of the ion In in the center-of-mass frame due to the 
potential v Va is 



3 

(x|^(+J;t) p - = -^C/e.L/^e^^*), (27) 



with an ion current density 

2 3 3 

fr" = jf-^s E E UejU^UZU* kj^-^-^-^t). (28) 
j,2=l ff=l 

If we sum over all flavors of the antincutrino and use the unitarity relation (I23|) . we obtain the 
total radial current density 

cr j=l 

As previously for two flavors, the modulations in (|28p disappear. 

However, the GSI experiments did observe modulations with an amplitude a w 0.18(3). Thus 
we must conclude that the unitarity relation f|23[) must be violated. Since so far only the lowest 
possible modulation frequency Alo = lj2~lo\ between the two lightest neutrinos has been measured, 
we may parametrize the right-hand side of the unitarity violation by 

KjUal = U 6jl + U 2 l(Sj,2Sl >1 + S jt x6i, 2 ) + • • • , (30) 



G 



and find 



3r = 



M„r 2 



{u S S,, 



+ 2u 2 iS 2 i cos[Akr - Aut + A0]} , 



(31) 



where 
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•21 = 



\U e2 U* el I (k 2 + fcO/2, A0 = arg U e2 U* el . 



(32) 



Assuming that the violation of unitarity is small, the sums Sq and S\ are close to unity. Then we 
deduce from tha experimental result a 0.18(3) that 



A possible origin of this unitarity violation could be that there are more than three families 
of leptons in nature and that universality of weak interaction is not valid for all of them. If the 
symmetry between quarks and leptons of the standard model persists to higher energies, we do 
not expect more than eight lepton families to exist — more than eight quark families would ruin 
asymptotic freedom and thus confinement. Thus there is room for more than the three quark and 
lepton families observed so far. Indeed, a fourth set of families is under intense discussion [l5| in 
connection with the new accelerator LHC at CERN. So far, there are only weak bounds on their 
masses from different sources H: 



If any of the heavier leptons is coupled with a coupling constant that does not fit into the CKM 
scheme, unitarity will certainly be violated. More data will be needed to decide precisely how. 



It is noteworthy that this analysis, in which we extract the properties of the unobserved an- 
tineutrino from the behavior of the ion, corresponds precisely to the usual entanglement analysis 
of decay processes such as 7r° — > 7 + 7. There the measurement of the polarization of one photon 
tells us immediately the polarization properties of the other, unobserved photon. 

A few comments are in place on several recent publications [lH E3, EH, EH, H3, HH which deny the 
relation between neutrino oscillations and the nonexponential decay seen in the GSI experiment for 
various reasons. In Ref. [lj|, the basic argument is that the antineutrino oscillations set in after 
their emission, so that they cannot be observed in the GSI experiment. The present discussion 
shows that although the forst part of this argument is true, the conclusion depends on the unitary 
assumption of the mixing matrix. 

Finally we should point out that similar oscillation phenomena in the associate production of 
particles together with an oscillating partner have been proposed and controversially discussed 
before by many authors in the production of muons together with antineutrinos in the decay 
7r~ — ► [i~ + [23j] , and in the production of A hyperons together with neutral Kaons [24|, In 
the latter case the oscillation would come from a nonunitarity of the quark mixing matrix, which 
seems to be much smaller than that of the neutrino mixing matrix reported here. 



Acknowledgment: We are grateful to Giinter Kaindl for organizing regular scientific evening talks 
in the Magnus-Haus Berlin of the German Physical Society. This brought us together on March 
11, 2008, in a lecture by P.K. on the GSI experiments, leading to the interpretation presented 
here. We also appreciate financial support by the WE-Heraeus-Stiftung for the 'Nachsitzung' in 
the Remise of the Magnus-Haus. 




(33) 



m t , > 256 GeV, m b , > 128 GeV, m T , > 100.8 GeV, m v , > 90.3 GeV. 



(34) 



VI. COMMENTS 
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VII. APPENDIX: PROPERTIES OF OUTGOING WAVE 

The interaction is time-dependent and we must adapt the scattering theory to this situation. 
Recall briefly the theory for a time-independent interaction, where the scattering amplitude is 
obtained from the standard limiting formula 

(p'|S|p) = lim (p'|p (+) (*))= lim (p'|Uj(*,0)|pW)= lim (p'|e^ ot e~ l( ^ 0+f)t |P (+) ) 

t— *00 t—*OC t—*OQ 

= lim e , '< JS p'- E »> t <p / |P (+) >- (35) 

t — >oc 

Here |p') denotes an eigenstate of the free Hamiltonian Hq with momentum p' and energy E p > = 
p /2 /2M, and |p^ + ^) is an eigenstate of the interacting Hamiltonian Hq + V with momentum p and 
energy E p . It solves the Lippmann-Schwinger equation: 

| P W) = |p> + 1 (36) 

E p - H + 17] 



which is verified by multiplying both sides by E — Hq from the left. Inserting this into (|35|) leads 
to 

(p'|5|p> = (p'|p) + lim f iEp '~ EP l)j p'\V\p^), (37) 
t— too — hi-pi + %r\ 

where r\ > is an infinitesimally number. The second term contains the T-matrix T p /. p = 
(p'|V r |p' + ') which describes true scattering. In the absence of neutrino oscillations, |p) is sim- 
ply the initial ion at rest, and (p'| the state with the ion Ijj with momentum p + k and the 
antineutrino v e with momentum — k. The limit t — > oo in the prefactor can simply be taken after 
rewriting it as —i dte lt - E p'^ Ep ^" > " ,t , which obviously tends to —2niS(E p > — E p ) in the limit 
t — > oo. The (^-function ensures the conservation of energy in the process. This is, of course, the 
standard derivation of Fermi 's Golden Rule which we repeated here to clarify that it is applicable 
only to processes in which the final state is an eigenstate of the free-particle Hamiltonian operator 

Another way of deriving this result is based on the spatial wave function associated with the 
state Ip^). One multiplies Eq. (f3"o| by the state (x| from the left and obtains the wave funcion 

(x|p<+>) = <x|p) + J dVG( J E p ;x,x')y(x')<x'|p(+)), (38) 



where 



1 C r1 3 n' pip' (*-*') pip' E r , ,. , 

G(E\ x, x') = (x| s |x') = / — 4t m -2M e ~ tp ^ , (39) 

V ' ' ' V l E-H +ir) J (27r) 3 J B-p' 2 /2M + iry Aur ' 1 ' 



with r = |x — x'|, x = x/r and p' E = \J2ME. 

In Born approximation, one inserts on the right-hand side of Eq. (|38p a plane wave (x'|p( + )) w 
(x'|p) = e lx ' p , and Eq. (J38]) becomes 

ip' E r p 

(x|pW) = (x|p) - 2M \ d 3 x'e-' i ( p '-P) x V(x'), (40) 

4nr J 

where p' is short for the momentum of the outgoing particle of energy E p in the direction of x: 
p' = pg p x. Thus, in Born approximation, the amplitude for the final particle to emerge with 
momentum p' is proportional the Fourier transform of the potential at the momentum transfer 
Ap = p' p. If the potential is a plane wave of momentum — k, i.e., if 



(41) 
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then the final state has the wave function 

pip'r 

(x|p(+)> = (x|p> - 2Mg ^ (3) (P' - p + k). (42) 

Let us adapt this formalism to the oscillating situation. According to Eqs. (fT3")) , (fT4"]) , the emission 
of an antineutrino of mass m\ and momentum — k is described by the time-dependent interaction 
potential 



e ikx 



v auc (x,t) = ^ M = ^k 2 + m 2 , (43) 

where we have dropped the factor cos 2 9 accompanying the coupling g, for brevity. As before, 
the incoming ion I has the momentum p. Its energy is p 2 /2M. The outgoing ion Ijj has the 
momentum p' and an energy E p > — Mh — M + p' 2 /2Mh- 

Let us first adapt the Lippmann-Schwinger approach. We introduce the time-dependent interact- 
ing state |p( +) (i)), which is an eigenstate of the full Hamiltonian, and satifies the time-dependent 
Schrodinger equation 

id t \p^(t)) = [H + v(x,t)}\P {+ Ht)), H = p 2 /2M. (44) 
The formal solution of this is 

|p (+) (*)) = U{t)\p), U^^fe-^o^'lHo+v^t)}^ ( 45 ) 
An implicit exression for this state can be written, by analogy with (|36|) . as 

|p<+>(t)) = |p) e -<* + 1 -V-(x,t)|pW(t)), E° p = p 2 /2M. ,-!(,! 

id t - Ho + vq 



This can again be verified by multiplication from the left with idt — Hq. Multiplying (|4l)|) by e lE p t , 
we obtain 

e^*|p (+) W> = IP) + -e^p*t>(x,t)e-^*|p(+)(0)), (47) 

To lowest approximation, we replace |p( + )(0)) by |p) and insert (|4"3")) to find 

/oo 
dt' G(t,t')e lk *e l ^ l{t '- to} \p)e-^- E >' , (48) 
-oo 

where G(t, t') is the Fourier representation of the operator (idt + E p — Hq + irj)^ 1 : 

r°° dE e -;#(*-*') 
G(t,t') = / ^ x (49) 

Performing the integral over t' in (|48| yields 

o , p-i(Bp-J5j-u> v ,i)* 

e iB p*|p (+) (£)) = |p) + g s e ikx |p)e- Mk ' lt0 . (50) 

E p - Ho- w k) i + if] 

We multiply this equation from the left and insert in front of the right-hand state |p) a completeness 
relation J d 3 a;|x')(x'| = 1. Approximating the matrix elements (x\(E p — Hq +i?7)~ 1 |x') as usual 
in the large-x regime by 

1 C d 3 n' P *p'(x-x') p ip' k r 



Ep - Ho - Wk,i + ir) J (2tt) 3 S p - p /2 /2M - w k ,i + i?7 ' 47rr 

where pj, is the momentum of the ion Ih which conserves the energy, i.e., p' k 2 /2M = E p — Wk.i- 
With this we obtain 

(xle^lp^ff)) = (x|p) - 2M J./y f d 3 a ./ e -i P ' fc 4x' e i(p+k)x' e -iu, k>1 t _ (52) 

(27r)' 3 47rr j 
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We now perform the integral over x'. For this we assume the inital state to have zero momentum, 
p + k = 0. The integral over x' forces the momentum of the outcoming ion Ih to be equal to k. 
The integral over x' creates a 5-function (2tt) 3 5^ (p' k x — p — k), so that we obtain 

(x|e ,;B p*|p(+)(t)) = (xlp) -2M-^-e-^ kx 'S^(p' -p + k)e-^ t0 . (53) 

4irr 

Note that since energy and momentum are balanced, then u k ,i — lo\ of Eq. ([5]). 

Consider now the case of two oscillating mass states and let us study the temporal behavior of 
the emerging energy distribution. The experiment does not explore the limit of very large times 
but measures the i-dependence starting from small t after the ion enters the storage ring. Instead of 
the limiting energy conservation 5-function —2TriS(E p i — E p ) in (|37[) . it observes an approximation 
to it valied for short times. To find it we insert, instead of (|43| . the mixed potential (fT4|) into 
Eq. (HT)) , so that the time-dependent factor in the resulting equation of type (j4"8)) has the form 



— i 



dt 



cos 2 9 e* E i.u-*W + sin 2 9 e *(^.v-^)* 



(54) 



where E^v = \/k' 2 + M' H + av,i and Eu — Mh + Q- Since m- < Q< Mh, we can approximate 
-Ei.k' " Ek « lu' — LOi where + tti 2 /2Mh- Let us write u)x,2 = |Am 2 /2M# = u> =p |Aw. 

Then (|5"4"|) becomes, with the abbreviations C = cos 2 and S = sin 2 0, 



/ „i(u)'—wi)t _ i i(w'-w 2 )t _ 1 

[C ; -+S-, 



M ' (w'-wi)/2 (w'-W2)/2 J 



(55) 



The absolute square of (|5"5|) mutiplied by some factor |T| 2 determines probability P(£) to find the 
initial ions in the ring at the time t. The integral over the final momenta is dominated by the 
immediate neighborhood of the poles at u>i t 2 where |k| = Q. There we may ignore the k-dependence 
of |T| 2 , approximating it by a constant, and obtain the probability for the ion Ih to emerge with 
an energy Mh + Q — u>' in the center-of-mass frame 

P u \t) sa {C 2 s 2 1 (uj') + S 2 s 2 2 (uj') + 2CScos{Aujt/2)s 1 (uj')s2{u')} \T\ 2 , (56) 

where Sj(u/) = sin[(cj' — u)i)t/2]/[(ui' — Wi)t/2]. For large t, the limiting relation sin 2 at/a 2 — » tTi8{a) 
allows us approximate s 2 (uj') ~ 2i:t5(u)' — cj,), and thus the first two terms in (|56[) by 

P£(t) pb 27ri[C* 2 (5(w' - wi) + S 2 8{u' - u 2 )]\T\ 2 . (57) 

If this is integrated over J d 3 k' / (2ir) 3 Q 2 J dui' /2ir 2 , the probabilities of v\- and ^-decays simply 
add, thereby yielding the ordinary /3-decay rate / — * Ih + v e without mixing. Consider now the 
third term in (|56|) . Here the integral over all u' yields 

J dw'Pg'(t) w 2C5cos(A^/2)2^^^^|T| 2 . (58) 

Thus we obtain for the total decay rate as a function of time 

P(t) = J dJ +/§''(*) ^2ir[l + 2CScos(Acut)]\T\ 2 . (59) 

It would be interesting to observe experimentally the predicted distribution (|56jl of antineutrino 
energies to' by measuring the recoil momenta k of the final ions Ih- The distribution consists of 
two peaks associate with the emission of the antineutrinos V\ and v-i- Centered between them lies 
the oscillating distribution proportional to Si(ui')s2(u)') shown in Fig. 1. 

Note that the usual Feynman diagrams in momentum space cannot be used to describe the 
observed oscillations as done in Ref. [17] . since they imply taking the limit t — > oo in which the 
oscillations disappear. Only diagrams in spacetime involving a propagator matrix in ^1,^2 field 
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Figure 3: Distribution of the temporally oscillating part of the antineutrino center-of-mass energies uj' in 
the decay I — > Ih + v e for different times t [the function si(u/)s2(u/) in Eq. ((56}]. 

space with off-diagonal matrix elements are applicable, and these reproduce the above-calculated 
oscillations. 
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